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Abstract. In this paper we describe a Fabry-Perot interferometer in the language of quantum optics. We go 
on to model the Fabry-Perot interferometer as a beam splitter having frequency dependent transmissivity 
and reflectivity coefficients. The antibunching, a totally non-classical effect, is to be expected if one excites 
this interferometer with carefully frequency-selected single photons arriving simultaneously from opposite 
directions. Contrary to a normal beam splitter, even slightly different frequency single photons should be 
able to show this effect, as long as the photon counters are not frequency selective. 

PACS. 42.50.-p Quantum optics 


1 Introduction 

Ever since its introduction, more than a century ago, the 
Fabry-Perot (FP) interferometer [1] has been ubiquitous 
in spectroscopy and high resolution interferometry. A cou¬ 
ple of decades later, laser resonators [2] adopted its work¬ 
ing principle. Other FP applications include gravitational 
wave detectors [3] and, more recently, for the needs of cav¬ 
ity quantum electrodynamics (CQED), ultra-high finesse 
superconducting [4] and fiber high finesse [5] Fabry-Perot 
cavities. 

The FP interferometer is an optical device composed 
of two highly reflecting mirrors and its classical descrip¬ 
tion [2] shows an optical device having extremely sharp 
variations going back an forth between total transmission 
and total reflection as the spacing between the two mir¬ 
rors (or the incident light frequency) is varied. There is a 
discrete set of frequencies, where, whatever the transmis¬ 
sivity of the individual mirrors is, the FP interferometer 
shows a total transmission. This effect can be explained 
classically, as well as quantum mechanically. 

The dynamic behavior of a FP interferometer has been 
investigated by Lawrence et al. in [6], with application 
in the context of Laser Gravitational-Wave Observatory 
(LIGO). It has been shown that a different behavior is 
observed in the case of a vibrational mirror and the ap¬ 
plication of a frequency modulated signal to the FP inter¬ 
ferometer. The extension of the model and experimental 
results of Lawrence et al. was done by Rohde et al. [7], 
where the modification of the natural decay time of atoms 
trapped inside the cavity was demonstrated. 
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Interesting experiments [8,9] showed that a single atom 
can behave like an optical FP cavity under certain con¬ 
ditions, thus bringing a traditionally macroscopic inter¬ 
ferometer to the atomic scale. In [10], Srivathsan et al. 
showed that a FP interferometer is able to reverse the 
exponentially-falling into an exponentially-rising tempo¬ 
ral envelope for a heralded single photon, a feature with 
potential interest in quantum information processing [11]. 
A recent proposal by Sun [12] uses FP interferometers to 
precisely resolve the continuous variable time-energy en¬ 
tanglement. 

Non-classical states of light have been used for decades 
and the most popular source of such states is the so-called 
spontaneous parametric down-conversion (SPDC) [13,14], 
where an incident pump photon, through a non-linear 
optical process, generates two output photons, generally 
called “signal” and “idler”. 

One non-classical effect is the so-called “HOM” or anti¬ 
bunching effect [15], experimentally demonstrated by Hong, 
On and Mandeb It manifests itself in a dip in the coin¬ 
cidence rate at the output of a beam splitter, when two 
single-photon Fock states impinge simultaneously at its 
inputs. 

The quantum optical description of a FP interferome¬ 
ter is typically done by quantizing the classical field modes 
[16]. In this paper we start directly from a quantum optical 
model of the FP interferometer. Being an optical lossless 
component, it is expected that it will show a SU{2) sym¬ 
metry [17,18]. 

Therefore, using a graphical method introduced in [19] 
and already used in [20] , we describe the input-output field 
operator transformations of a FP interferometer. Next, we 
go on to predict an antibunching effect, showing a more 
rich structure, due mainly to the heavy frequency depen- 
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dence of the transmissivity and reflectivity coefficients of 
a FP interferometer. 

This paper is organized as follows. In Section 2 the 
quantum optical description of the FP interferometer is 
sketched. Its equivalence to a beam splitter having frequency- 
dependent transmissivity and reflectivity coefficients is em¬ 
phasized, too. The prediction of the antibunching effect 
when two single-photon Fock states impinge on the Fabry- 
Perot interferometer from opposite directions is discussed 
in detail in Section 3. The extension to the case when the 
input light is non-monochromatic is done in Section 4. The 
conclusions drawn in Section 5 close this paper. 


2 The quantum optical description of a 
Fabry-Perot interferometer 

A Fabry-Perot interferometer (depicted in Fig. I) is com¬ 
posed of two parallel (sometimes slightly concave) highly 
reflecting mirrors placed at a distance L. Throughout this 
paper, we assume the mirrors to be identical^ and we re¬ 
strict our analysis to a single axis (i.e. orthogonal to the 
mirrors). Moreover, we assume each highly reflecting mir¬ 
ror to be frequency independent in the optical range of 
interest. Thus, we can model it by a beam splitter (BS) 
having a transmissivity T and a reflectivity R. Since we 
assume our mirrors to be of negligible thickness compared 
to the distance L, we take the transmissivity to be a pos¬ 
itive real number^ i.e. T = e and the reflectivity to be 
purely imaginary, therefore, we have 

T = e and R = i\/ 1 — (I) 

One can easily check that these coefficients obey the well- 
known relations for a beam splitter jTp -|- |i?p = 1 and 
TR* + T*R = 0 [21]. The angle of incidence of our light 
beam of interest is denoted by 9 (see Fig. 1). 

The FP interferometer is resonant at wavelengths Aq 
where the round trip in the cavity of length L causes a 
phase shift Lp = kL that is an integer multiple of 27r. 
Therefore, for 6 —>■ 7r/2 we have the relation L = NXq/2 = 
ttIVc/wo, where ujq is the resonant (angular) frequency of 
the cavity, TV € N and c is the speed of light in vacuum. 
Since k = uj/c one gets the phase shift ip = ttNuj/ujo = utg 
where we denoted to = Ljc. This expression will be used 
later on. 

Since the interesting case from a practical point of view 
is when the angle of incidence is almost normal {9 —> 7r/2), 
we shall consider this scenario for the rest of the paper. 
In this case, our FP interferometer can be considered as 
having only four modes (ports) as depicted in Fig. 2 (see 
also the discussion in Appendix A). We can find the input- 
output operator transformations by various means. In Ap¬ 
pendix A, they are obtained via a graphical method. The 

^ The case with non-identical mirrors was discussed in [16]. 
Extension to this case is also possible with the graphical 
method discussed in Appendix A. 

^ For high reflectivity FP cavity mirrors, e is a small real 
number. However, in all computations that follow, this condi¬ 
tion is not needed. 



Fig. 1. The Fabry-Perot interferometer. Each input/output 
mode is labelled with a number, starting with 0. We associate 
to each input/output mode I a field creation operator denoted 
by aj. 



Fig. 2. The creation field operators for a Fabry-Perot interfer¬ 
ometer when 0 —>■ 7r/2. Input modes are labelled by I and 
V, while output modes are labelled by T and R. Each in¬ 
put/output mode has the corresponding creation operator at¬ 
tached to it. 


end results are 

a\ = Tfphlp -\- ( 2 ) 

and 

Qy ~ RfpO^j' (3) 

where Tfp and, respectively, Rfp are given by equations 
(30) and, respectively, (31). 

By analyzing equations (2) and (3) one notes that these 
are actually the operator transformations of a beam split¬ 
ter having frequency dependent transmissivity 

rp2„itjjto 

Tfp (w) = ^ _ jp2^i2ojto 


and reflectivity 


X _j_ gj2ajto 

Tfp (w) = ^ ^2gi2wto 


( 5 ) 


coefficients. The probability of transmission and, respec¬ 
tively, reflection of a monochromatic beam of light of fre¬ 
quency u! is easily computed as 


\Tfp{co)\^ 


1 

\ + cos^ (wto) 


( 6 ) 
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Fig. 3. The Fabry-Perot interferometer can be also seen as a 
beam splitter with frequency-dependent Tfp (uj) and Rfp (u) 
(transmission and reflection coefficients). 


and 


lo /,.m2 T^COS^Ujto) 


( 7 ) 


where we replaced T and R according to equation (1) 
and the cavity finesse was defined as = 2\/l — je^. 
We note that we have |T/p(a;)p -b |i?/p(a;)p = 1 and 
Tfp (uj) Rfp {u}) + T^p (w) Rfp (w) = 0 i.e. we have a SU{2) 
symmetry, characteristic of linear lossless optical systems 
[17,18], 

Therefore, the whole Fabry-Perot interferometer can 
be modelled as a simple beam splitter (see Fig. 3), with 
its transmissivity and reflectivity coefficients being heavily 
frequency dependent. 

It is then expected that, similar to a beam splitter, 
the same classical and non-classical effects from quantum 
optics must be reproduced by a FP interferometer. 


3 Antibunching with a Fabry-Perot 
interferometer in the case with 
monochromatic photons 

The antibunching is a non-classical effect observed when 
two single-photon Fock states impinge simultaneously on 
the two inputs of a balanced (50/50) beam splitter [15]. 

Since the Fabry-Perot interferometer is formally equiv¬ 
alent to a beam splitter (with frequency dependent trans¬ 
missivity and reflectivity, as discussed before), using the 
jl/ly) input state should show this nonclassical effect at 
its output. Therefore, we apply the input state to a FP 
interferometer 

\tpin) = [l/ly) = 4«y|0) (8) 

where both input light quanta are supposed to be mono¬ 
chromatic of frequency to. Using again the input-output 
field operator transformation relations (2) and (3) takes 
us to the output state vector 

IV'out) = |0) 

= TfpRfp {\ 2 t 0 r) + lOT2fi)) + {Tfp + R)p) IItIa) (9) 

We shall denote by chom (e, w) = Tfp (uj) + Rfp (w) and 
Prom (Ejw) = \chom {£,^)\'^ the amplitude, and, respec¬ 
tively, the probability of coincidence detection at the out¬ 
puts T and R of the FP interferometer. 



transmittivity of the individual mirrors, s 


Fig. 4. The antibunching factor Phom = \chom\^ versus the 
transmissivity of individual mirrors e and the phase variation 
LutoiTT. The blue thick curve corresponds to pairs (e, lj) yielding 
PHOM{e,ijj) < 0 . 01 . 


In order to obtain the antibunching effect, we impose 
the amplitude of the llrlfi) state to be zero. The condi¬ 
tion Chom (w) = 0 implies cos^ (wto) = s'^/4, (l — £^), an 
equation that can be readily solved yielding 

1 f \ 

UJ+/- = —arccos ±— , (10) 

to \ 2VT^J ^ ’ 

This equation has solutions as long as e < Eq, where eg = 
2 {y/2 — l). We could see eg as the maximum value of 
the transmissivity of the individual mirrors, so that the FP 
interferometer can still exhibit the antibunching effect^. 

In practice, however, having Phom (e, w) below a cer¬ 
tain value would be still satisfactory. In Fig. 4 we plot 
Phom{£,>-^) in respect with the transmissivity of individ¬ 
ual mirrors (e) and the phase variation {ujtg /tt) . The thick 
blue curve from Fig. 4 outlines the values e and uj for 
Phom < 0 . 01 . 

Up until now, these results are equivalent to the ones 
obtained with a balanced beam splitter. The true interest 
in using a FP interferometer lies in the fact that the same 
effect can be obtained using photons of different colors. 

We consider now two (monochromatic) photons of two 
different frequencies, Wg and Therefore, the input state 
can be written as 

|'0m)=^ da;dw'C/(w)Cy(w')al(w)a|/(a;')|0) (11) 

with Cl {ui) = S {u! — ujs) and Cy(w) = S{uj — uii). The 
frequency-dependent input-output creation operators now 
yield 

C'j (^s) — Pfp (^s) ^fP (^s) 

® It is interesting to note that at the frequency given by 
equation (10), we have Tfp (a;+) = 1/2 (1 -|- i) and Rfp (a;+) = 
l/2( —l-|-i) while at the frequency uj- we have Tfp{ijj-) = 
— 1/2 - I - i/2 and Rfp (ui-) — 1/2 -|- i/2. 

















4 


S. Ataman: The QO description of a Fabry-Perot interferometer and the prediction of an antibunching effect 


and 

(^0 ~ ^fp (^0 d” (^0 (d^) 

The output state vector is obtained by considering the 
field operator transformations (12—13) on \ tjjin) from equa¬ 
tion (11) and gives 

\' 4 ^out) — ^fp (^s) ^fp (^z) (^t) | 0 ) “b 

^fp ^fp (^s) (^s) |0) “b 

(^Tjp (cl’s) ^fp (^i) (^s) i.^i) 

“b Rfp iy^s) Rfp ij^i) (^0 (^s)^ |0) (Id) 

Since the photon counters are typically frequency non- 
selective (at least in a small spectral range), a coincidence 
detection at the T and R outputs can be modelled as 

CHOM (ws,Wi) = (0| JJ dujduj'ax {uj)dR (w')IV’oixt) 

= Tfp (ws) Tfp (uji) + Rfp (ws) Rfp (wi) (15) 
where we used the commutation relations [21] 

[oj (w) ,aj'(w')] = ^ (w - w') (16) 



Fig. 5. The pairs of frequencies (uJsjUJi) satisfying 
From (uJajtOi) < 0.01. The thick red curve corresponds to 
e = 0.7, the green one corresponds to e = 0.4 and the blue 
curve corresponds to £ = 0.1. 


where I G {T,R}. We shall be interested in the second 
order correlation function [21,24] 


for I G {T,R}. The antibunching condition now reads 
chom (wsjWi) = 0, leading to the constraint on Wg and 

OJi'. 


(t, t + x) = {i/jout\ e!^ (t) E^j^\t + r) 

E^+'^ {t + T)E^+^ {t)\ilJout) (20) 


cos {ujsto) cos {ojito) = — -(17) 

4 (1 - 

Again, in a practical scenario, we might relax our con¬ 
straints and, instead of imposing a zero amplitude for 
Chom (wgjWi), we could ask for a value oiPnoM (^ 5 ,^^) = 
\chom (ws, Wi) P below a certain limit. The points (ujs^uji) 
satisfying Phom (wg,Wi) < 0.01 are plotted in Fig. 5 for 
three different values of e. 

In Appendix B we discuss what changes implies the 
more realistic assumption that the signal and idler pho¬ 
tons originate from a SPDC process. 


where e\~^ {t) = (<)]^ and t (t+r) denotes the detec¬ 

tion time at the output T (R). We will use the Schrodinger 
picture therefore we shall evolve JV’m) to \ipout)- Using 
again the field operator transformations (2) and (3) ob¬ 
taining the output state vector 

lipout) = jj diodj (Pfp (w) (w) -b Rfp (w) a|j (w)^ 

^ (j^fp {oJ ) a^x {w ) + Tfp (w ) a|j (w )^ C, (w, w ) jO) (21) 

Combining now equations (19), (20) and (21), after some 
computations, one ends up with^ 


4 Antibunching with a Fabry-Perot 
interferometer in the case of 
non-monochromatic photons 


G‘TM + r) = T 


dwdu;'e-“‘e-*‘^'(‘+^) 


2 


We extend now the previous analysis to non-monochromatic 
photon wave packets [22,23]. The input state can be ex¬ 
pressed as 

\ipin) = jj dwdw'C(w,u;')a7(w)av(w')|0) (18) 


^ (^/p (^) '^fp ) T Rfp (w) Rfp {lo )) C (w, to ) 


( 22 ) 


where we used again the commutation relations (16). If 
we make the simplifying assumption that the signal and 
idler frequencies are not frequency correlated i.e. we have 
C (w, w') = ^0 (w) Cl equation (22) simplifies to 


where C (w, ui') is the bi-photon wave packet and it is nor¬ 
malized via // 1C pdwdo;' = 1. The output positive 

frequency electric field operators have to be extended to 

(19) 


G^^^ (t, 1 -b r) — jC(f {t) Cf (t -b t) -I- C(f {t + t) Cf (t) 1^(23) 

^ Typically, for the antibunching effect we shall be interested 
in (t) which is the time integrated version of G^^^ {t, t + x) 
i.e. G^^^ (r) = f G^^^ {t, t + x) dt. 


eI+^ it) = ^l dG2e-“‘az (oi) 
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where we have the inverse Fourier transforms 

Ql (<) = ^ I dt^e—‘T/p (..) C™ (ca) (24) 

and 

c (i) = ^ I dt^e—(w) Cm (C^) (25) 

with m G {0,1}. To the author’s best knowledge, there is 
no closed-form solution to these Fourier integrals. How¬ 
ever, one can easily check that if Co (w) —>■ 5 (w — Ws) and 
Cl (w') —(5 (w' — uJi) equation (22) migrates into 

^ \Tfp (lOs) Tfp (uJi) + Rfp (lOs) Rfp (wi) | (26) 

If uJs and oji obey equation (17), then the antibunch¬ 
ing effect is assured. If one assumes the SPDC-like fre¬ 
quency correlation, C (w, w') —?► d (w — Ws) <5 (w' — uip + uig) 
the discussion from Appendix B still holds. 


5 Conclusions 

In this paper we started from the full quantum optical 
description of a Fabry-Perot interferometer and showed 
that it can be modelled through a beam splitter having 
heavily frequency-dependent transmission and reflection 
coefficients. 

Owing to this model, the antibunching effect, well- 
known for beam splitters has been predicted for the Fabry- 
Perot interferometer, too. However, even different frequency 
photons can show this effect, this being in contrast with 
the normal antibunching. Monochromatic, as well as non- 
monochromatic input light was considered. 


A The computation of the field operator 
transformations in a FP interferometer 

In the following we will employ the graphical method in¬ 
troduced in [19]. Beam splitters will be depicted by the 
butterfly-like structure and delays will simply add a e*”^ 
factor. The graphical representation of the Fabry-Perot in¬ 
terferometer is detailed in Fig. 6. Each crossing of a highly- 
reflecting mirror is modelled by a beam splitter (having 
the transmission and reflection coefficients T and, respec¬ 
tively, R) and each passage inside the cavity yields an ex¬ 
tra factor of where ip = kL and k is the wavenumber 
of the (monochromatic) light that illuminates the inter¬ 
ferometer. 

It is not difficult to find the expression of the input 
creation operator a\ in respect with the output ones. We 
can differentiate between transmission type creation oper¬ 
ators (dg, djy ...) and reflection type operators (d^, 
dg, a\/^ ...). Adding up all these contributions, one ends 
up with 

+R (dj -b -b ...) (27) 


As discussed before, the interesting case from a practical 
point of view is when the angle 6 tends to 7r/2. The many 
inputs and outputs tend to group together in the following 
way: the inputs (d|, dy, djg ...) become indistinguishable® 
and we relabel them as a single input /, whose creation 
operator will be denoted by a\. Similarly the reflected out¬ 
puts {a\, dg, a\^ .. .) become also indistinguishable and 
will be denoted by d|j. The same logic applies in regroup¬ 
ing the appropriate transmission operators (dg, dg, djg, 
djy ...) into dy. The “dark” ports (dj, a\, djg, djg ...) 
are grouped into the “vacuum” port dy (see Fig. 2). Now 
equation (27) takes the much simpler form 


d| = T^e*"^ (1 -b -b -b ...) d^ 

-bi? (1 -b -b + ■ • ■) 4 (28) 


Since we have now two trivial geometric series and |i?| < 1, 
equation (28) can be written in the rather simple form 



TZTJ^^t + R 



T2^2v \ 

1 - i72e*2v ) 


(29) 


We could define the transmission 

fP 1 _ Jl2^^2<p 


(30) 


and reflection 


Rfp = R 



T2Qi2ip \ 

1 - i?2e*2vp j 


(31) 


coefficients of the FP interferometer allowing the easy 
writing of the field operator transformation given by equa¬ 
tion (2). Similar arguments allow one to write dy in re¬ 
spect with dJp and dj^, as done in equation (3). 

The regrouping of input and output field operators 
that became indistinguishable allows a much simpler graph¬ 
ical representation of the FP interferometer, as depicted 
in Fig. 7. 

The infinite series of beam splitters from Fig. 6 was 
replaced by a loop, having the nodes B, C, D and E in 
Fig. 7. The meaning of the loop is that besides the direct 
path, there is a possibility to loop ones, twice and so forth. 
For example, in Fig. 7, the direct path from d^ to d\ took 
simply a factor of e*"^ (besides the factor caused by 
elements outside of the loop). Looping ones implies, be¬ 
sides e*''’, a supplementary factor of . Looping twice 

implies a supplementary factor of and so forth. 

This is the geometric series (1 -b R^e^^^ + R'^e'^'^P -b ...) 
we found before and it trivially yields 1/(1 — R'^e^^P). In 
the end, from the recursive graphical method we obtain 
T^e®‘^/(1 — R^e^'^P), a result identical to equation (30). 

Therefore, as a rule, in the graph from Fig. 7, every 
input-output path having a common segment with this 
loop must take a supplementary factor of 1/(1 — R^e^'^P). 

® Although one might wish to input a photon only at the 
port 1, the extreme closeness of the ports 1, 7, 13 ... makes 
impossible to tell which port actually received the photon. 
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T T 


Fig. 6. The graphical description of a Fabry-Perot interferometer for light incident at an arbitrary angle. Each passage of the 
beam through the highly reflecting mirrors is modelled with a beam splitter. The labelling was done in accordance to Fig. 1. 




Fig. 7. The graphical description of a Fabry-Perot interferome¬ 
ter in the recursive form when light impinges at a normal angle. 
The blue circles denote the input ports while the red ones de¬ 
note the output ports. Every path having a common segment 
with the loop B — C — D — E gains a factor of 1/(1 — 


For example, the input creation operator d| can be 
reached from the output port d|j through the reflection 
coefficient of the input beam splitter yielding a factor R 
(and touching no loop). There is a second contribution 
from d|j to a\: transmission T from the first (left) beam 
splitter, a factor of e*"^, a R (from the second i.e. right 
beam splitter), another factor of e®"^ and finally a factor 
of T to d|. This path has common segments with the loop, 
therefore it also takes the factor 1/(1 — Adding 

up the two possible routes gives R + /{I — R'^e^'^'^) 

and we found straight away the result from equation (31). 


B The antibunching condition with photons 
originating in a SPDC process 


In the SPDC process, the signal and idler frequencies are 
bound by the relation ojp = ujs + tOi where ojp denotes the 
pump frequency. This constraint modifies equation (17) 
and takes us to a single variable (say Wg) yielding 


ap cos {2oJsto) + Y® 1 - sin {2uJsto) = Pe - ctp (32) 

where we denoted ap = cos (ujpto) and = e^/2 (l — e^). 
Using well-known trigonometric identities, we transform 
the sines and cosines into tan(u;sto) and have 


/3etan^ (wgto) - (wgto) + Pe - ‘2-ap = 0 (33) 

This second order equation in tan(a;sto) yields two solu¬ 
tions 


tan(a;gto) \+/- 


al± al 

JT 


/3| -b 2apl3e 
-(34) 


as long as the quantity under the square root is not neg¬ 
ative. This constraint imposes 

al - 2ap/3e -b - 1 < 0 (35) 

This can happen only if ap is 

Pe-l<ap<P, + l (36) 

The second inequality is always satisfied since /?£ > 0, 
therefore we are left with 

- 1 < cos (wpto) (37) 

This relation puts constraints on the values of to, uip and 
e in order to make possible the antibunching effect at the 
signal frequency tUg. 
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